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Abstract 

Let y be a vector space over a field or skew field F, and let U 
be its subspace. We study the canonical form problem for bilinear or 
sesquilinear forms 

U {V/U) X 1/ ^ F 

and linear mappings U ^V,V ^U,V/U V, V ^ V/U. We solve 
it over F = C and reduce it over all F to the canonical form problem 
for ordinary linear mappings W ^ W and bilinear or sesquilinear 
forms W X W ^ ¥. Moreover, we give an algorithm that realizes 
this reduction. The algorithm uses only unitary transformations if 
F = C, which improves its numerical stability. For linear mapping 
this algorithm can be derived from the algorithm by L. A. Nazarova, 
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1 Introduction 

In this paper, we give canonical matrices of bilinear or sesquilinear forms 



[/ X 1/ ^ C, 



{V/U) X 1/ ^ C, 



where \^ is a complex vector space and U is its subspace. 

We use the following canonical matrices of bilinear or sesquilinear forms 
on a complex vector space given in pQ (see also |3] ) . Two square complex 
matrices A and B are said to be congruent or *congruent if there is a nonsin- 
gular S such that S^AS = B or, respectively, S*AS = B, where S* := S'^ 
denotes the complex conjugate transpose of 5*. Define the n-hj-n matrices 







-1 



1 1 
-1 -1 
1 









1 

1 i 



A 1 
A 





Theorem 1 ([H p. 351]). (a) Every square complex matrix is congruent to a 
direct sum, determined uniquely up to permutation of summands, of matrices 
of the form 

J„(A) o" 

-1)"+-'^, and A is determined up to replacement by A^^. 



in which A 7^ 0, A 7^ 
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(b) Every square complex matrix is *congruent to a direct sum, determined 
uniquely up to permutation of summands, of matrices of the form 











in which |A| = 1 and 
matrix A„ instead of r„,. 



> 1. Alternatively, one may use the symmetric 

□ 



A canonical form of a square matrix for congruence/*congruence over any 
field F of characteristic different from 2 was given in [5j up to classification 
of Hermitian forms over finite extensions of F. 

Let us formulate the main result. For generality, we will consider matrices 
over any field or skew field F with involution a ^ a, that is, a bijection on 
F such that 

01 + (3, 



a 



a 



a + /5 = a + /3, a(3 = Pa, 
for all a,f3 e¥. 

We denote the m-by-n zero matrix by Omn, or by Om iim = n. It is agreed 
that there exists exactly one matrix of size n x and there exists exactly one 
matrix of size x n for every nonnegative integer n; they represent the linear 
mappings — F" and F" — and are considered as the zero matrices 0„o 
and Oon- For every p x q matrix Mpq we have 



pq 







mO 



'M„ " 






OpO 






. 0^0. 






OmO_ 







and 



pq 



0, 



On. 



M 







pq 

Oo„ 



M^pq Opfj 



%q 



On 



[Mpq %n] . 



In particular. 

For each matrix A 







'pO 







0, 



pq- 



\a 



over F, we define its conjugate transpose 
A* 



A^ 



\a 



If S*AS = B for some nonsingular matrix 5*, then A and B are said to 
be ^congruent (or congruent if F is a field and the involution on F is the 
identity — in what follows we consider congruence as a special case of *con- 
gruence). 
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A sesquilinear form on right vector spaces U and V over F is a map 

g-.u ^¥ 

satisfying 

Q{ua + u(3, v) = aQ{u, v) + j3Q{u , f ), 
Q {u, va + v'f3) = Giu,v)a + G (m, v')(3 

for all u, u' G ?7, f , v' G and a, /3 G F. If F is a field and the involution 
on F is the identity, then a sesquilinear form becomes bilinear — we consider 
bilinear forms special case of sesquilinear forms. 

If Cl, . . . , Cm and fi,...,fn are bases of U and V, then 

Gef = [aij], aij := G{ei,fj), (1) 

is the matrix of Q in these bases. Its matrix in other bases e'^, . . . , and 
f [,..., f^ can be found by the formula 

Ge'f = S*GefR, (2) 

where 5* and R are the change of basis matrices. 
For every u & U and v & V, 

g{u,v) = [u]lGef [v]f, 

where [u]e and [v]f are the coordinate column- vectors of u and v. 
In this paper, we study sesquilinear forms 

U xV (V/U) X 1/ ^ F, (3) 

in which ?7 is a subspace of so we always consider their matrices ([1]) in 
those bases of U and V that are concordant as follows. 

Definition 2. Let Q be one of sesquilinear forms ([3]), in which V is a. right 
space over F, and U is its subspace. Choose a basis ei, . . . ,en of V such that 

ei, . . . , Cm is a basis of U if Q : U x V ^ ¥, 
Cm+i, . . . , e„ is a basis of U if Q: {V/U) xV—*¥. 
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By the matrix of Q in the basis Ci, . . . , e„, we mean the block matrix 



[A\B] 



ami 



ttl,m+l 



a 



m,m+l 



air 



(5) 



in which 



a,; 



^(e„e,) if^:f/xy^F, 



'Ai 





Bi 


" 





A2 





B2 



By the block-direct sum of block matrices and [yl2|-B2]) we mean 

the block matrix 



In Section O we will prove the following theorem (a stronger statement 
was proved in [21 Theorem 1] in the case U = V). 

Theorem 3. Let ¥ be a field or skew field with involution {possibly, the 
identity if F is a field), V be a right vector space over ¥, and U be its 
subspace. Let Q be one of sesquilinear forms 



UxV 



(V/U) xV ^¥. 



(6) 



(a) There exists a basis Ci, ...,€„ of V satisfying (jl]), in which the matrix 
of Q is a block-direct sum of a p-by-p matrix 



[K\Opo], K is nonsingular, 
and matrices of the form 

[JMM (Ol), [JM\E,] (g^o), 



in which 






1 



En := 0, 



01 



(7) 



(8) 



(9) 
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{the summands ^ or ^ may he absent). The block K is determined by Q 
uniquely up to *congruence, and the summands of the form are determined 
by Q uniquely up to permutation. 

(b) // F = C, then one can replace in this direct sum the summand ([7]) 

by 

where © ■ • ■ © is the canonical form of K defined in TheoremUland each 
Ki is Pi-by-pi. The obtained block-direct sum is determined by Q uniquely up 
to permutation of summands, and so it is a canonical matrix of the sesquilin- 
ear [in particular, bilinear) form Q . 

Let us formulate an analogous statement for matrices of linear mappings. 

Definition 4. Let F be a field or skew field, be a right vector space over 
F, and U be its subspace. Let A be one of linear mappings 



V/U ^V, V ^ V/U. 



Choose a basis Ci 



, . . . , e„ of such that 



ei, 



e„ is a basis of U, 



ifU-^VoiV^U, 



(10) 



Cm+i, . . . , e„ is a basis of U, if V/U ^ V or V ^ V/U. 

By the matrix of A in the basis ei, . . . , e„, we mean its matrix in the bases 

'ei,...,e^ of [/, ifU^VoiV-^U, 
ei + U,...,em + U of V/U, if V/U ^ V oi V ^ V/U, 



and ei, 



, Cn of V. We divide Ag into two blocks 
A 



A, 




if t/ y or V/U V, 
iiV -^U oiV V/U, 



(11) 



where A is m-hj-m. 

The following theorem will be proved in Section [51 
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Theorem 5. Let ¥ be a field or skew field, V be a right vector space overW, 
and U be its subspace. Let A be one of linear mappings 

U^V, V^U, V/U^V, V^V/U. (12) 

(a) There exists a basis ei, . . . , e„ of V satisfying flTU]) . in which for the 
matrix Ae of A we have: 

A^, ifU^VorV/U^V, 
Ae, ifV^UorV^V/U 

is a block-direct sum of a p-by-p matrix 

[K\Opo], K is nonsingular, (13) 

and matrices of the form 

[J,(0)|0,o], [Umi (14) 

where Eg was defined in ([9]) {the summands ( fT3l) or (fT4l) may be absent). The 
block K is determined by A uniquely up to similarity, and the summands of 
the form (1141) are determined by A uniquely up to permutation. 

(b) If ¥ = C, then one can replace the summand f|T3|) by a block- direct 
sum of square matrices of the form 

The obtained matrix is determined by A uniquely up to permutation of sum- 
mands, and so it is a canonical matrix of the linear mapping A. 

We do not rate Theorem [5] as new; it is readily available from the canon- 
ical form problem solved in [U § 2] . We include it in our paper since the 
singular indecomposable summands of the canonical forms in Theorems [3] 
and [5] coincide, and our proofs of Theorems [3] and [5] are similar and are based 
on regularization algorithms that decompose the matrix of each form ([3]) and 
each mapping ( |T2l) into a block-direct sum of 

• its regular part [K\Opo] with nonsingular K (see ([7]) and (fT3|) ). which is 
determined by ([3]) or (fT2!) up to *congruence or similarity, and of 

• its singular summands of the form [Jg(0)|Ogo] and [J(,(0)|-E'q] (see ([8]) 
and dH])), which are determined uniquely. 
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If F = C, then these algorithms can use only unitary transformations, which 
improves their numerical stability. These algorithms extend the regular- 
ization algorithm [2] for a bilinear/sesquilinear form, which decomposes its 
matrix into a direct sum of a nonsingular matrix and several singular Jordan 
blocks. An analogous regularization algorithm was given by Van Dooren [7] 
for matrix pencils and was extended to matrices of cycles of linear mappings 
in 0. 

The canonical form problems for matrices of forms ([3]) and mappings (fT2!) 
are special cases of the canonical form problem for block matrices, whose form 
resembles 



© 




Definition 6. By a bangle over F we mean a matrix 

A=[A^\...\ Ak-i I Ak I Ak+i I ... I A] 



(15) 



over F, partitioned into vertical strips, among which one strip Ak is square 
and boxed. The number of rows of A and the number Ui of columns of 
each strip Ai are nonnegative integers. Let 



B = [Bi I ... I Bk-i 



Bk 



B, 



k+l 



...IB,] 



(16) 



be another bangle with the same sizes of strips and the same k and t. We say 
that the bangles A and B are *congruent or, respectively, similar and write 



A^ B or A^ B 
if there exists a nonsingular upper block-triangular matrix 

. . . Sit 



(17) 



S 







Stt 



{Sii is Hi X Hi) 



over F such that 



B 



Sl^AS or B = Sj;i^AS. 
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Then 

Bk = Slp,AkSkk or Bk = Sf^^AkSkk, 

this means that the boxed strips of *congruent/similar bangles are *congru- 
ent /similar. The following lemma is obvious. 

Lemma 7. Two bangles are *congruent/ similar if and only if one reduces to 
the other by a sequence of the following transformations: 

(a) Any transformation with rows of the whole matrix, and then the 
^congruent/ similar transformation with columns of the boxed strip {this 
transformation reduces ( 1151) to 



[EA^\... I EAkE* I EAu+i | . . . | EAt] 

or, respectively. 



[EAi\... \EAk-i I EAkE-^ \ EA^+i | . . . | EAt] 
with a nonsingular E) . 

(b) Any transformation with columns of an unboxed strip. 

(c) Addition of a linear combination of columns of the i^^ strip to a column 
of the j^^ strip if i < j . □ 

Note that the canonical form problem for matrices of forms ([3]) and map- 
pings (fT2l) is the canonical form problem for bangles (fTSi) with two strips. But 
applying our algorithm to bangles with two strips we can produce bangles 
with three strips (see Section 13.21) ; so we consider bangles with an arbitrary 
number of strips. 

The paper is organized as follows. In Section [2] we formulate our main 
theorem about the existence of a regularizing decomposition of a bangle. In 
Sections [3] and H] we construct regularizing decompositions of bangles with 
respect to *congruence and similarity. In Section [5] we use these decomposi- 
tions to prove the main theorem and Theorems [3] and [51 
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2 Bangles 



In this section, we formulate our main theorem, which reduces the canonical 
form problem for bangles up to * congruence/similarity to the canonical form 
problem for nonsingular matrices up to * congruence/ similarity, and solves it 
for complex bangles. 

By the block-direct sum of two bangles f|T5|) and f|T6l) with the same num- 
ber of strips and the same position of the boxed strip, we mean the bangle 



A\£B :-- 



Ai 
Bi 



Ak 











At 
Bt 



Definition 8. A regularizing decomposition of a bangle 
A= [Ai\...\Ak- 



Ak 



A 



k+l 



over a field or skew field F with respect to *congruence/similarity is a bangle 
Syi satisfying two conditions: 

(i) T,A is *congruent/similar to A, and 

(ii) is the block-direct sum of 

— its regular part 



[Opo I . . . I Opo K Opo I • • • I Opo] , K is nonsingular, (18) 



and its singular part being a block-direct sum of matrices of the 
form 



[0,o|...|Ogo ^,(0) 0,o|...|0,o] 



[■■■\E,\... J,(0) ...], [... J,(0) 



(19) 
], (20) 



in which Eg is defined in ([9]) and the dots denote sequences of 
strips Oqo- 

Both the regular and the singular parts may have size 0-by-O. 
The following theorem generalizes Theorems [3] and [5l 
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Theorem 9. (a) Over a field or skew field ¥, any bangle A possesses regu- 
larizing decompositions for "^congruence and for similarity, their regular parts 
are determined by A uniquely up to ^congruence and, respectively, similarity, 
and their singular parts are determined by A uniquely up to permutation of 
summands. 

(b) If ¥ = C and is a regularizing decomposition of a bangle A for to 
*congruence, then its regular part flTS]) is ^congruent to the block- direct sum 



1+) [Op,o I ... I Op,o [k^ I ... I Ko] 



in which Ki ® ■ ■ ■ ® Kg is the canonical form of K defined in Theorem [I] 
and each Ki is Pi-by-pi. Replacing in the regular part by this block-direct 
sum, we obtain a canonical form of A for ^congruence {in particular, for 
congruence) since the obtained bangle is ^congruent to A and is determined 
by A uniquely up to permutation of summands. 

(c) // F = C and is a regularizing decomposition of a bangle A for 
similarity, then its regular part is similar to a block-direct sum of matrices 
of the form 



[0,o|...|Ogo ^g(A) 0,0 |... I 0,0], A^O. 



Replacing in the regular part by this block-direct sum, we obtain a canon- 
ical form of A for similarity since the obtained bangle is similar to A and is 
determined by A uniquely up to permutation of summands. 

Note that for bangles with respect to similarity this theorem can be de- 
duced from the canonical form problem solved in [4, §2]. 



3 Regularization for ^congruence 

We give an algorithm that for every bangle over a field or skew field F con- 
structs its regularizing decomposition for *congruence. If F = C, then we can 
improve the numerical stability of this algorithm using only unitary transfor- 
mations. The algorithm is the alternating sequence of left-hand and right- 
hand reductions, which we define in Sections 13.11 and 13.21 
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3.1 Left-hand reduction for *congruence 



Let 



A=[A^\...\A, 



k~l 



Ak 



A 



k+l 



(21) 



be a bangle over F. Producing *congruence transformations (a)-(c) from 
Lemma [7] with A, we can reduce its submatrix [y4i|742 | . . . | ^fc-i] by the 
following transformations: 

(a') arbitrary transformations of rows; 

(b') arbitrary transformations of columns within any vertical strip Ai] 

(c') addition of a linear combination of columns of the i^^ strip to a column 
of the j^^ strip if i < j. 



(22) 



using transformations (b') with Ai and (a'), then make zero A22, ■ ■ ■ , ^2,A:-i 
by transformations (c'). Transforming analogously the submatrix 



First we reduce [Ai ^2 . . . 


Ak-1 


to the form 




" 


A12 




Ai^k-i 




/ 


A22 




A2,k-1 



A12 1 . . . I 



k-1 



we reduce fl22l) to the form 















B3 











/ 








I 












and so on. Repeat this process until obtain 



Bk-i 





" 
















" 






























^rk-1 






























r2 > 0, 



Tk ^ 0, 



(23) 



and extend the obtained partition into horizontal strips to the whole bangle 
( I2T1) . Make zero all horizontal strips of the blocks A^, . . . ^At except for the 
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first strip and obtain 



Ck{M) :-- 








/. 





Ir. 







Ml 


M2 . 


. Mk 





0^2 


. 





. 


• 0., 



Mfc+1 






Mt 






(24) 



(we have divided the boxed block Aj. into k vertical strips conformally to its 
partition into horizontal strips) for some 



M 



Ml 



Mo 



Mt] =: L{A). 



(25) 



Clearly, r2, . . . , are uniquely determined by A. 



Definition 10. We say that a bangle A reduces to a bangle B by admissible 
permutations and write 

A^ B 

if A reduces to i? by a sequence of the following transformations: 

• permutation of rows of the whole matrix and then the same permuta- 
tion of columns of the boxed strip, 

• permutation of columns in an unboxed strip. 
Clearly, 



A^B 



(in the notation (fT7|l ). 

Lemma 11. (a) The equivalence 

Ck{M) ~ A(iV) 

holds for all 

M = [\ Ml |M2|...|Mt], 



Ar^ B and A ~ 5 



M ^ N 



N 



iVi 



No 



and each k ^t. 

(b) // F = C, then for every bangle A we can find 
transformations. 



(26) 

using only unitary 
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Proof, (a) The equivalence fl26l) is trivial if /c = 1. Let k ^ 2. Reasoning by 
induction on /c, we assume that 



Ck-i{M) ~ Ck-i{N) ^ M^N 
and prove the equivalence fl2B]) as follows. 
(^) Suppose Ck{M) ~ Ck{N), that is, 

5^,/:fc(M)5 = A(iv) 

for some nonsingular 

Sn ■ ■ ■ Sit 
S= -.. : 

Since both Ck{M) and Ck{N) have the same first vertical strip 





(27) 



(28) 



(29) 



L 



(we join its zero horizontal strips), by fl28l) we have 



kk 






" 


Sii — 


'0 


" 













and so Skk has the form 



Si 



kk 



Let 



R 



Pi P2 
Ps 



S22 ■ ■ ■ S; 



(30) 



2f 



Su_ 

be a submatrix of ( 129|) with S'^fc of the form ( l30l) . Due to ( l28l) . 

P*A_i(M)/? = A_i(iV). (31) 
So /:fc_i(M) ~ /:fe_i(A^), and by ([27D M ~ A^. 
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=) Suppose M ^ N. By ([27]), C,k-i{.M) ~ Ck-i{N), this ensures 
for some nonsingular 



P 



11 



It 



Ptt 

Denote by Bi and Ci the strips of Ck~i{M) and £fc_i(A^): 
A_i(M) = [fill... |fife_i 



Bk+1 
Ck+i I 



C.]. 



Then 

A(M) : 

and by (I3l 
A(M) ^ 












Bk-i 






Bk+1 









Bk+2 









-P/i;+l,A;+l 



A.(M) 



■/ 

(-ffc+i.fe+i) 






p 



" 





Ci 




Ct ' 









C[ 









Bt 




where CJ, 



C'^ are some matrices. 



This proves ( l26l) . Let us give an alternative proof of ( l26i) using *congru- 
ence transformations (a)-(c) from Lemma [71 Due to that lemma, it suffices 
to show that those transformations (a)-(c) with ([2^ that preserve all of its 
blocks except for Mi, . . . , Mj produce all transformations (a)-(c) with (!25l) . 

• We can add a column of Mi to a column of M, if i < j. Indeed, in the 
case j ^ k this is a column-transformation within the boxed block of 
Ck{M), and so we must produce the *congruent row-transformation — 
add the corresponding row of the i^^ horizontal strip of f[Ml) to the row 
of the j^^ horizontal strip. This spoils zero blocks of the j^^ horizontal 
strip, but they are repaired by additions of columns of Ir ■ 
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We can also make arbitrary elementary transformations with columns 
of Mi if 2 7^ 1: in the case i ^ k these transformations spoil J^. but it 
is restored by transformations with its columns. 



(b) Let F = C. We must prove that if 
A=[A^\...\ Ak-i 



A. 



A 



k+l 



is reduced to (125|) by the algorithm from this section, then r2, . . . ,rfc and 
Mi,...,Mi can be found using only unitary transformations with A. By 
unitary column-transformations within vertical strips Ai,...,Af,_i of A 
and by unitary row-transformations, we sequentially reduce its submatrix 
[y4i|y42 I ... I Afe-i] to the form 



" 
















" 



































* 


* 















* 


* 







* 


* 




* 


* 



where each H^^ is a nonsingular rj-by-rj block and all *'s are unspecified 
blocks (this reduction was studied thoroughly in [6l Section 4]). The matrix 
A takes the form 








H, 





H, 



r2 



Ml 


M2 . 




* 




* 


* 


* 


■ *rk 



Mk+i 



Mt 



(32) 



in which *r 



■ ■ ,*rk ^2 y. r2, . . . ,rk X matrices. Replacing Hr^, 
by the identity matrices of the same sizes and all *'s by the zero matrices, 
we obtain (!25|) because (132|) can be reduced to (p5|) by those transformations 
(a)-(c) from Lemma [7] that preserve r2, . . . , and Mi, . . . , Mt. □ 



3.2 Right-hand reduction for *congruence 



Let 



A= [\~Ar]A2\...\At] 



(33) 
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be a bangle over a field or skew field F. 

First we reduce A by * congruence transformations 



[ SAiS* SA2 I ... I SAt ] , S* is nonsingular, 



to the form 










^3 




Bt+i 




B[ 


B'2 


B's 




B't+i 



(34) 
(35) 



in which the rows of [B[ B'^ are linearly independent and B2 is square. 
Then we make zero B'^, . . . , B[j^^ adding columns of B[ and B2, and as in 
sequentially reduce [_B3|S4 | . . . | Bt+i ] to the form 



" 















































Irt-1 













In 

















obtaining a partition of the first horizontal strip of (!35|) into t substrips. 
Conformally divide the first vertical strip of the boxed block into t substrips 
and obtain 



7^(M) 



On 


. 


































Or, . 




























Ir2 





.. 


• o,,_, 

















Irt~i 













.. 





On 








In 
















Ml 


M2 .. 


. Mt-i 


Mt 


Mt+i 























(36) 



for some 

M = [Mi\...\Mt 

with Mf+i = B2. 

Lemma 12. (a) The equivalence 

n{M) ~ 7^(iV) 



Mi 



t+i 



] =■■ R{A) 



(37) 



M ^ N 



(38) 
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holds for all 



M= [Ml 



Mi 



t+i 



], N=[Ni 



Ni 



t+i 



]■ 



(b) If ¥ = C, then for every bangle A of the form fl33l) we can find fl36l) 
using only unitary transformations. 

Proof, (a) Let us prove the equivalence (l38l) using *congruence transforma- 
tions (a)-(c) from Lemma [7] (alternatively, one could use induction on t as in 
the proof of Lemma [TlTa)). Due to Lemma [71 it suffices to show that those 
transformations (a)-(c) with ( l36l) that preserve all of its blocks except for 
Ml, . . . , Mt+i produce all transformations (a)-(c) with ( 137|) . 

• We can add a column of Mj to a column of Mj if i < j; by the definition 
of *congruence transformations we must add the corresponding row of 
the i*^ horizontal strip of ( l36l) to the row of the j*^ horizontal strip; 
although this spoils a zero block of the j*^ horizontal strip if i 7^ 1, but 
it can be repaired by additions of columns of J^^ . 

• We can also make arbitrary elementary transformations with columns 
of Mj if i ^ t: these transformations spoil J,., if i 7^ 1, but it is restored 
by transformations with its columns. 



(b) Let F = C. First we reduce the bangle (!33|) by transformations ([34 
with unitary S to the form (!35|) . in which the rows of [B[ B2] are linearly 
independent and B2 is square. 

Then we sequentially reduce [i?3|i?4 | . . . | -Bt+i ] by unitary column- 
transformations within vertical strips and by unitary row-transformations 
to the form 



' 




































Hr2 











Hrt-1 




* 


* 







* 


* 




* 


* 



where each Hr^ is a nonsingular rj-by-rj block and the *'s are unspecified 
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blocks. The matrix A takes the form 



Or, 


. 


. 








0^2 


. 








. 


• 0,, 





Ml 


M2 . 


. Mt 










H, 



n 





Hr. 



(39) 



where M^+i = i?2- Replacing iJ^j, . . . ,Hr^, by the identity matrices of the 
same sizes and all *'s by the zero matrices, we obtain (136|) because (139!) 
can be reduced to (!36|) by those transformations (a)-(c) from Lemma [7] that 
preserve ri, . . . , and Mi, . . . , M^+i. □ 

3.3 Regular izat ion algorithm for *congruence 

For any bangle 

A= [Ai|...|Afc„i| Ak \Ak+i\...\At] (40) 



over F, its regularizing decomposition for *congruence can be constructed as 
follows. 

Alternating the left-hand and the right-hand reductions for *congruence, 
we construct the sequence of bangles 

A':=L{A), A":=R{A'), A'" := L{A"), A"" := R{A'"), . . . 

until obtain 



= [^Opo|...|Opo] or = [Opo|...|OpoE] 



(41) 



with a nonsingular K. 

Producing this reduction, we in each step have deleted the reduced parts 
of A; say, in step 1 we reduced A to the form fl2^ and took only its unreduced 
part A' = L{A). Let us repeat the reduction of (140|) preserving all the reduced 
parts of A: 



In step 1 we transform A to CkiA') of the form fl24|) . 

In step 2 we reduce its subbangle A' to TZ{A") preserving the other 
blocks of Ck{A'), and so on. 
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After n steps, instead of fl41l) we obtain some bangle A, which is *congruent 
to A. Due to the next theorem, A is a regularizing decomposition of A up 
to admissible permutations of rows and columns. 

Theorem 13. If A is a bangle over afield or skew field ¥, then A reduces by 
admissible permutations of rows and columns to a regularizing decomposition 
of A for ^congruence. 

Proof. We give a constructive proof of this theorem. 

By admissible permutations of rows and columns, A reduces to a block- 
direct sum of the bangle f|T8|) in which K is the same as in fHTl) . and a bangle 
D in which each row and each column contains at most one 1 and its other 
entries are zero. We obtain a regularizing decomposition of A for *congruence 
replacing D in this block-direct sum by from the follows statement. 

Let D be a bangle in which each row and each column 

contains at most one 1 and the other entries are zero. 

Then D reduces by admissible permutations of rows and (42) 

columns to a block-direct sum of bangles of the form 

(UnD and (EOj). 

Let us prove (142|) . By admissible permutations of rows and columns of D, 
we reduce its boxed strip to a direct sum of singular Jordan blocks. Then 
we rearrange columns in each unboxed strip such that if its j) and 
entries are 1 and i < i' , then j < j'. It is easy to see that the obtained bangle 
So is a block-direct sum of bangles of the form f|T9|) and fl20|) : each singular 
Jordan block Jp(0) in the decomposition of Dk gives the summand f|T9l) if 
those row of D that contains the last (zero) row of Jp{0) is zero, and the 
summand fl2U]) otherwise. The summands (I2U1) with p = give zero columns 
in unboxed strips oi D. □ 

4 Regularization for similarity 

We give an algorithm that for every bangle over a field or skew field F con- 
structs its regularizing decomposition for similarity. If F = C, then we can 
improve the numerical stability of this algorithm using only unitary trans- 
formations. 
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4.1 Left-hand reduction for similarity 



Let 



A=[A^\...\A, 



k-l 



A 



k+l 



be a bangle over F. Using similarity transformations with A, we can reduce 
its submatrix | • • • | ^fc-i] by transformations (a')-(c') from Section 

I3.1[ We reduce this submatrix to the form 



' 



























Ir2 


















































ri > 0, 



Tk-l > 0, 



and obtain a partition of the bangle A into k horizontal strips. Then we 
divide the boxed block Ak into k vertical substrips of the same sizes, make 



zero all horizontal strips in the blocks Ak, 
and obtain 



, At except for the last strip, 



A(M) 



Ir, 















On . 








. 


■ o.,_, 





Ml . 


. M,,„i 







Mfc+1 





Mt 



for some 



M 



M, 



[Mi|...|Mfc_i 
Lemma 14. (a) The equivalence 

Ck{M) ~ Ck{N) 

holds for all 

M = [Mi\. 
[Ni\. 



Mk+i\...\Mt] =: LiA). 



M ^ N 



N 



|Mfc_i 
Nk-i 



Nk 



Nk+i I . 



...\Mt]. 
..\Nt]. 



(43) 



(44) 



(b) If ¥ = C, then for every bangle A we can find ( l43l) using only unitary 
transformations. 



21 



Proof, (a) This statement follows from Lemma [7] since those transformations 
(a)-(c) with that preserve all of its blocks except for Mi, . . . ,Mt produce 
all transformations (a)-(c) with (jH]). For example, we can add a column of 
Mi to a column of Mj if i < j: although in the case j ^ k we must subtract 
the corresponding row of the j^^ horizontal strip of P3l) from the row of the 
^th iiorizontal strip, and this may spoil zero blocks of the i^^ horizontal strip, 
but they are repaired by additions of columns of 7^.^. 

(b) Let F = C. By unitary column-transformations within vertical strips 
of A and by unitary row-transformations, we sequentially reduce its subma- 
trix [74i|y42 I • • • I ^fe-i ] to the form 



" 




* 


* 




* 


* 











Hr2 




* 


* 










































where each is a nonsingular rj-by-r^ block and all *'s are unspecified 
blocks. The matrix A takes the form 



Hr, 






in which * 



H, 




rk-i 







* 


* 


* 


*rk-l 


* 


Ml . 


• Mk-l 





Mfc+1 



* 



(45) 



5 *rk- 



are ri x ri,...,rfc-i x Vk-i matrices. Replacing 
by the identity matrices of the same sizes and all *'s by the 
zero matrices, we obtain (HSl) since (HSil reduces to (H3|l by those transforma- 
tions (a)-(c) from Lemma [7] that preserve ri, . . . , Vk-i, Mi, . . . , Mf. □ 



4.2 Right-hand reduction for similarity 



Let 



A= [\Ar\A2\...\At] 

be a bangle over F. 

First we reduce A by similarity transformations 



(46) 



[ SAiS SA2 I ... I SAt ] , S" is nonsingular. 



(47) 
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to the form 







B2 


i?3 




Bt+i 










B', 




B't+i 



(48) 



in which the rows of [-Bi i?2] are hnearly independent and Bi is square. 

Then we make zero -B3, . . . , Bt+i adding columns of Bi and i?2, and 
sequentially reduce [ i?3 | . . . | B[j^^ ] to the form 



" 















" 






































Irt 























The matrix A transforms to 
7^(M) = 



Ml 


M2 


M3 . 


. Mt 


Mt+i 

















Or, 


. 


. 






















0,3 . 


. 














Jj.3 








. 


• 0,, 























. 


. 

















for some 



M 



Ml M2 I . . . I M^+i ] =: R{A) 



with Mi = Bi. 

Lemma 15. (a) The equivalence 

n{M) ~ 7^(iV) 

holds for all 

-iMi+i], 



M ^ N 



M =[ Ml 



Mo 



N =[ Ni 



No 









Ir, 





(49) 



(50) 



(b) // F = C, then for every bangle A of the form ( H6l) we can find ( l49i) 
using only unitary transformations. 
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Proof, (a) It is easy to show that those transformations (a)-(c) from Lemma 
[7] with (1491) that preserve all of its blocks except for Mi, . . . , M^+i produce 
all transformations (a)-(c) with fl50l) . Say, we can add a column of Mi to a 
column of Mj if i < j: although we must subtract the corresponding row of 
the j^^ horizontal strip of fH9l) from the row of the i^^ horizontal strip, and 
this spoils zero blocks of the i^^ horizontal strip if j 7^ t + 1, but they are 
repaired by additions of columns of J^^. 

(b) Let F = C. First we reduce A by transformations (147|) with unitary S 
to the form (HHj) . in which the rows of [Bi B2] are linearly independent and 
Bi is square. 

Then we sequentially reduce its submatrix [S3IS4 | . . . | iJ^^^ ] by uni- 
tary column-transformations within vertical strips and by unitary row- 
transformations to the form 



' 


Hr2 


* 


* 




* 


* 















* 


* 



















Hrt 























where each Hr^ is a nonsingular rj-by-rj matrix. The matrix A takes the form 



Ml 


M2 . 


. Mt 







Or, . 


. 








. 


• 0,, 





0., 


. 


. 





K 



r2 








Hr 




(51) 



by the identity matrices of the same sizes and all *'s by 



Replacing Hr^ , 

the zero matrices, we obtain fH^ since fl^ reduces to fH^ by those transfor- 
mations (a)-(c) from Lemma [7] that preserve r2, . . . , r^+i. Mi, . . . , M^+i. □ 



4.3 Regularization algorithm for similarity 

For any bangle 

A = [Ai\...\Ak-i \ Ak \Ak+i\...\At] 



over F, its regularizing decomposition for similarity can be constructed as 
follows. 
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First we apply subsequently the left-hand reduction for similarity to A 
until obtain 



Bk 



5, 



fc+i 



L{L...{L{A))...)= [0^o|...|0^o 
in which the first k — 1 strips have no columns. 

Then we apply subsequently the right-hand reduction for similarity to 

Bk+i I . . . I -Bi ] 



B 



Bu 



until obtain 



Rs{Rs...{Rs{B))...)= \\ K |0„o|...|Ono1 
with a nonsingular K. 



(52) 



Producing this reduction, we in each step have deleted the reduced parts 
of A. Let us repeat the reduction preserving all the reduced parts of A 
and denote the obtained bangle by A. Clearly, A is similar to A. Due to 
the next theorem, A is a regularizing decomposition of A up to admissible 
permutations of rows and columns. 

Theorem 16. If A is a bangle over a field or skew field F, then A reduces by 
admissible permutations of rows and columns to a regularizing decomposition 
of A for similarity. 

Proof. We give a constructive proof of this theorem. By admissible permuta- 
tions of rows and columns, A is reduced to a block-direct sum of the bangle 
f ITS]) with K from flS^ and a bangle D in which each row and each column 
contains at most one 1 and the other entries are zero. Replacing D in this 
block-direct sum hj T^d from P2|) . we obtain a regularizing decomposition 
of A for similarity. □ 



5 Proofs of Theorems M, S], and [5 

Proof of Theorem (a) Let us prove the statement (a) for *congruence; its 
proof for similarity is analogous. 

Let A be a bangle over F. In view of Theorem [131 A possesses a regular- 
izing decomposition for *congruence, which is obtained from A by admissible 
permutations of rows and columns. 
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Let El and S2 be two regularizing decompositions of A. Then Ei ~ S2. 
We need to prove that 

Sf^~E^2'^ and Sf^^Ef^ (53) 

where and are the regular and the singular parts of Sj (i = 1, 2). 
If 

L(Si) = i?(Si) = El, 

then El = Si*"^ and ([53]) holds. 

Let Iv(Si) 7^ Si or -R(Si) 7^ Si. Suppose for definiteness that 

i^(Si) ^ Si. (54) 

Each row and each column of S^™^ {i = 1, 2) contains at most one 1, the 
other entries are zero. Due to this property, the reduction of Sj to 

a := Ck{L{J:,)). (55) 

of the form fl24p can be realized by admissible permutations: 

Si ^ fi.; (56) 

moreover, L(Sj) is a block-direct sum of a bangle of the form f|T8|) and a 
bangle, in which each row and each column contains at most one 1, the other 
entries are zero. By fH2]) . L(Sj) reduces by admissible permutations of rows 
and columns to its regularizing decomposition, so we may take Qi such that 
L(Sj) is a regularizing decomposition. 

Since Si ~ S2, we have i^i ~ ^2, and so by (l26l) and (|55|) 

L(Si) ~L(S2). 

Due to flMl) . the size of Iv(Si) is less than the size of Si, reasoning by induc- 
tion we may assume that flS5]) holds for -L(Sj); that is, 

L(Si)''^s ~ L(S2)'''s and L(Si)^'°s ^ L(S2)''''^. 

Then 

ni'^^ni'^ and fifs^(]fs 
since ^i and ^2 have the form fl24l) . This proves fl53l) due to fl56|) . 

(b) This statement follows from (a) and Theorem [TJ 

(c) This statement follows from (a) and the uniqueness of the Jordan 
Canonical Form. □ 
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Proof of Theorem O Let Q be one of sesquilinear forms 

UxV (V/U) xV -^¥. 

Let us prove that the canonical form problem for its matrix [A\B] (defined 
in (l5l)) is the canonical form problem under *congruence for the bangle 



[[A]b] or [b[J]], 

respectively, and so Theorem [3] follows from Theorem [91 

It suffices to prove that a change of the basis of V reduces [A\B] by 
transformations 



[A B] 



^ < 



S*[A B] 



S*[A B] 



S P 

Q 

s o' 

P Q 



lig-.u xv 



F, 



(57) 



if ^: {V/U) X y ^ F, 



in which 5* and Q are nonsingular matrices and P is arbitrary. 
Case 1: [A\B] is the matrix of 

g-.uxv^w, ucv, 

in a basis Ci , . . . , e„ of satisfying (jl]) . If 

fj = ^iPli ^ dnpnj, j = l,...,n, 



(58) 



is another basis of V such that /i, . . . , is a basis of U, then the change 
matrix from Ci, . . . , e„ to /i, . . . , /„ has the form 



R = [P^, 



'S P 
Q 



where S is the change matrix from ei, . . . , Cm to /i, 
the matrix [A\B] reduces by transformations (ISTj) . 

Case 2: [A\B] is the matrix of 



fm in U. Due to ©, 



^: (V/t/) X V 



F, 



f/ C K 
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in a basis ei, . . . , of satisfying (jl]). If (135]) is another basis of V such 
that fm+i, • • • , /n is a basis of U, then the change matrix from ei, . . . , e„ to 
fi,...,fn has the form 

-R = [Pij] ~ p Q 

where S is the change matrix from ei + U, . . . ,em + U to fi + U, 
in V/U. Hence, the matrix [A\B] reduces by transformations ([5 

Proof of Theorem El Let A be one of linear mappings 

U^V, V-^U, v/u-^v, v^v/u. 

Let us prove that the canonical form problem for its matrix 

A 



,fm + U 
□ 



A. 



B 



if t/ ^ r or V/U 

[A\B] liV oiV ^V/U, 
see (ITTj) ) is the canonical form problem under similarity for the bangle 



B' 



A' 



A' 



B' 



or 



respectively, and so Theorem [5] follows from Theorem [9l 

It suffices to prove that a change of the basis of V reduces Ae by trans- 
formations 



'A 








'A 


B 




. Q~\ 




B 



V, 



[A B] 
[A B] 



S-^[A B] 
S-^[A B] 



S * 

Q 

S 

* Q 



'A 




's-^ " 




'A 


B 




* Q'\ 




B 



S ifA:U 

iiA:V- 

ifA:V^ V/U, 
S HA: V/U ^V, 



(59) 
(60) 
(61) 
(62) 



in which S and Q are nonsingular matrices and the *'s denote arbitrary 
matrices. 
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Case 1: is the matrix of 



A-.U-^V or A:V^U, U C V, 
in a basis ei, of V satisfying (ITUi) . If 

fj = (^iPij ^ enPnj, j = l,...,n, 



(63) 



is another basis of V such that fi,...,fm is a basis of f/, then the change 
matrix from ei, . . . , e„ to /i, . . . , has the form 



R = [Pi 



'S P 
Q 



where 5* is the change matrix from ci, . . . , to /i, 
matrix A^. reduces by transformations fl5^ or (IBUl) . 

Case ^; is the matrix of 



fm in t/. So the 



A:V/U^V or ^: V ^ V/f/, f/ C V, 



in a basis ei, . . . , e„ of \^ satisfying (UOl) . If ( ]63l) is another basis of ^ such 
that /m+i, • ■ • , /n is a basis of f/, then the change matrix from ei, . . . , e„ to 
fi,...,fn has the form 

r 1 [5 0" 



R 



P Q\ ' 

where S is the change matrix from ei + U, . . . ,em + U to /i + f/, ...,/„ + f/ 
in Hence, the matrix A^ reduces by transformations fl6T|) or fl62l) . □ 
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